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The running of Newton's constant can be taken into account by considering covariant, non local 
generalizations of the field equations of general relativity. These generalizations involve nonanalytic 
functions of the d'Alembertian, as ( — D)~ a , with a a non integer number, and ln[— □]. In this paper 
we define these non local operators in terms of the usual two point function of a massive field. We 
analyze some of their properties, and present specific calculations in flat and Robertson Walker 
spacetimes. 
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I. INTRODUCTION 

In quantum field theory, the coupling constants depend on the energy scale according to the Rcnormal- 
ization Group Equations (RGE). This running behavior is crucial to understand pcrturbative and nonper- 
turbative aspects of elementary particles and their interactions. 

Over the years, there have been several attempts to incorporate the running behavior in astrophysics and 
cosmology. In the first case, a scale dependent Newton's constant would explain at least part of the dark 
matter. There are proposals in the literature [l[ where it is assumed that the gravitational constant around 
a compact object depends on the radial coordinate r. The dependence is dictated by the RGE, replacing 
the energy scale ft by 1/r. 

In a cosmological framework, the expansion of the universe makes any physical length scale to increase 
with time and naively one would expect the coupling constants to become time dependent as the universe 
evolves. Hopefully, this time dependence would help to explain the accelerated expansion of the Universe and 
to shed light into the dark energy problem. This point of view has been considered by a number of authors, 
and several proposals have been analyzed in detail. For example, it was assumed [H that Newton's and 
cosmological constants depend on t according to the RGE through the replacement ft — > 1/i. The Einstein 
equations with these time dependent constants in general do not satisfy Bianchi identites, and therefore 
some constraints are needed for mathematical consistency [3|. Related proposals assume that the correct 
replacement is obtained with the Hubble "constant" jl — > H [j] , the inverse radii of the event horizon [5[ , or 
with the scale factor p, — > ji/a(t) Q. Ressumation of the Schwinger DeWitt expansion in the large curvature 
limit produce running couplings that depend on the scalar curvature 0- 

In field theory, the running of the coupling constants is computed in momentum space, and the relevant 
momentum scale is associated to the characteristic energy in a scattering process. However, in astrophysics 
and cosmology, one is eventually interested in the running of the couplings in configuration space. As a 
trivial property of Fourier transforms, a scale dependent coupling constant in momentum state implies in 
general a non local dependence in configuration space. Therefore, if the running of couplings is produced by 
quantum effects (they could also be produced by extra dimensions, see below), in principle one should be able 
to compute a non local effective action and obtain real and causal effective equations for the gravitational 
field. 

This program can be explicitly worked out for weak gravitational fields, in the case where quantum fields 
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are free and massless. Indeed, after integrating the quantum degrees of freedom, the effective action [8| is 
equivalent to the classical action with the replacement of the coupling constants by the non local kernels Q 



ai(n) -> «;(-□) = a i0 {pt) + Qln 



□ 



(1) 



where fi is a reference energy scale and the values of the constants Cj are fixed by the RGE. Therefore, in this 
particular example the running of the gravitational constants is correctly and completely taken into account 
by replacing p 2 = — □ at the level of the effective action (this approach has also been considered in Refs. 
(lOl[ll|). The effective field equations for the gravitational field give rise to a quantum corrected Newtonian 
potential, that in turn can be interpreted as the classical one with a running gravitational constant. 

As already stressed, if one assumes that the modifications to General Relativity come from quantum 
effects, the running behavior of the coupling constants should be contained in the effective action [l2|], along 
with other quantum corrections. The "Wilsonian" approach, in any of the versions described above, is a way 
of partially taking into account the effects of quantum fields. The running of the coupling constants is not 
necessarily logarithmic. For example, renormalization group analysis based on the effective average action 
of quantum Einstein gravity [l3| and non perturbative studies in the lattice theory of quantum gravity [l4| , 
suggests different behaviors for the scale dependence of the gravitational constants, as for instance a power 
law scaling of Newton's constant. One could implement this running in either cosmological or astrophysical 
scenarios following the ideas described above, by an adequate replacement of the momentum by r _1 or t . 
A covariant approach that encompasses both scenarios, based on the use of the kernel £?(—□) has also been 
suggested [15| . This is difficult to implement at the level of the Einstein-Hilbert action, because formally 
any (analytic) function of the dAlembertian acting on the Ricci scalar is a total derivative and does not 
contribute to the field equations. When implemented at the level of the equations of motion, the modified 
Einstein equations take the form 



- nRg^> = 8wG{-u)T f 



fXV 



(2) 



For example, in the context of lattice theory of quantum gravity, it has been suggested that the strength 
of gravitational interactions might slowly increase with distance as a consequence of vacuum polarization 
effects. Provided that a non trivial ultraviolet fixed point exists, to leading order in the vicinity of such fixed 
point the scale dependence of the Newton's constant is characterized by a critical exponent v = (2a) -1 , and 
the size of the corrections is set by a non perturbative scale L. In this scenario, the manifestly covariant non 
local gravitational coupling can be written as fl5| 



G(-D) = G N (1 + a Q L- Aa (-U + L ) 



(3) 



where Gjv is the usual Newton's constant, clq is a positive number whose typical value was estimated to 
be of order of 10 -2 , and estimates for the value of v^ 1 = 2a vary from v^ 1 w 3.0 to v~ x « 1.7 (l4j |. The 
length scale L separates the ultraviolet regime where non perturbative corrections can be neglected from 
long distance regime where such corrections become significant. This scale might be related to the vacuum 
expectation value of the curvature of the spacetime, and therefore^Jn cosmological situations, it is generally 
chosen to have a typical value of order of the Hubble ratio H~ 



15] 



Similar proposals have been considered to introduce infrared corrections to General Relativity, which may 
help to understand the present acceleration of the universe. These are in turn inspired by extra dimensional 
brane world models (see below) [l6j . 

Effective equations like Eq.([2]) are in general inconsistent, i.e. the Bianchi identity is not satisfied. There 
is a way out, which consists in considering the effective field equations in the weak field limit. In this case, 
it is possible to write a non local effective action whose variation gives the effective field equations to linear 
order in the curvature [l7j . This procedure is consistent, although restricted to weak fields. A more drastic 
approach is to assume that Eq.([2]) is valid to all orders, and that the allowed energy-momentum tensors are 
only those that satisfy (G(-n)T^ v ). u = 0. 
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Another problem of Eq.© is that, unless the non local corrections contained in G{— □) are treated 
perturbatively, the theory will contain ghosts. Indeed, as shown in Ref . [l8j , at the linearized level the 
ghost-free theories must be of the form 

G$ ~ m 2 (-D)(h^ - VlUf h) = 8irG N T^ (4) 

where G^jj is the linearized Einsten's tensor, g^ v = r]^ v + and h = rj^h^. Eq.Q generalizes the Pauli- 
Fierz theory of massive gravitons to the case of a non local mass m 2 (— □). Note that the linearization 
of Eq.([2]) is not of the Pauli-Ficrz form. 

In brane world scenarios, the non local kernels previously considered in the literature involve cither arbi- 
trary powers or the logarithm of the d'Alcmbcrtian. For example, in order to modify General Relativity at 
distances larger that a certain scale L, it has been assumed that the leading behavior of the kernels is 

Gn m 2 {-U) 



G(-u) 



L- 2a (-D)- a (5) 



where a is a positive number [2(| (a = 1/2 for the five dimensional DGP model [lj|) and L is a length scale 
which for cosmological situations is naturally assumed to be of order of H^ 1 . Note that in contrast to the 
running suggested by non pcrturbative analysis of quantum effects (see Eq. [3]), the Newton's constant here 
decrease with distance. As an illustration, a "phcnomcnological" running coupling of the form 

with bL <C 1, interpolates between the usual Newton's constant at laboratory distance scales and G(0) <C Gn 
at scales much larger than L. 

In order to analyze the effective field equations it is necessary to define properly the action of an arbitrary 
power of the d'Alembertian on scalar and tensor functions. This is one of the main goals of the present 
paper. We will introduce a formal definition based on an integral representation, as it has been previously 
done with the logarithm of the d'Alembertian Q, and use this representation to obtain the action of the 
non local operator on arbitrary functions. We will consider both astrophysical and cosmological situations. 
This problem has been previously considered by Hamber and Williams, who computed (—□)"" by analytic 
continuation of the action of (—□)", for an integer number n > 1, to the case n — —a € R, in Robertson 
Walker [l5[ and Schwarzschild metrics [2l|. However, this approach has some drawbacks. On the one hand, 
it is only possible to get an answer for the action of (— n) _Q on some particular functions that depend either 
on t or r. On the other hand, the non local and causal properties of the kernel are not apparent from its 
definition. 

The paper is organized as follows. In the next Section we introduce our formal definition for (— n)~". 
In Section [HII we calculate it explicitly in flat spacetime, and reobtain previous results based on a different 
definition. We also show that the distribution ln[— □] can be obtained as a particular limiting case. In Section 
IIVI wc compute explicitly the action of (— Q on a scalar function independent of the spatial coordinates 
in Friedman Robertson Walker (FRW) spacetimes. We consider the de Sitter spacetime, and other FRW 
metrics with initial singularities. We also study the action of the kernel ( — □ + b 2 )~ a on a constant function, 
and check whether it gives b~ 2a or not, a result that has been used to argue that Eqs. like © can explain 
the smallness of the cosmological constant. In Section|V]we use the results obtained for the kernel ( — D)~ a 
to derive similar ones for ln[— □]. As a concrete application, in Section IVT1 wc discuss the effect of a running 
Newton's constant on the generation of gravitational waves. We summarize our results in Section [VIII 

Throughout the paper we set fl = c = 1 and adopt the sign convention denoted (+++) by Misncr, Thornc, 
and Wheeler 123 . 
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II. FORMAL DEFINITION OF (-□)" 



We will define the kernel (— □) a using an integral representation in terms of the two point function for 
massive fields. The definition is based on the formal identity 



(-□)-" = 



2 sin(a7r) 



+00 



dmm^ 2a - 



1 



- □ + m 



valid for < a < 1. A similar representation has been previously considered [8|] to define ln[— ^-], i.e. 



In 



□ 
J 2 



4-00 



dm 



1 



-□ + m 2 



(7) 



(8) 



) Q , it is possible to obtain ln[- 



It is worth noting that both kernels arc related: knowing the kernel (- 

by considering the identity x~ a = 1 — a In x + 0(a 2 ) (Tlj . 

The tensorial character of the kernels is defined by the one of the two point function of the massive field. 
When (— D)~ Q acts on a scalar (tensorial) function, (— □ + m 2 ) -1 represents the propagator of a scalar 
(tensorial) massive field. In what follows we will mainly consider the action on scalar functions. 

The definitions given above are unambiguous in Euclidean space, where — □ + m 2 is a unique positive 
definite operator. In Lorentzian spacetimes, it is necessary to impose an additional prescription to select 



one of the propagators (— □ 



') (retarded, advanced, Feynman, etc). As we will be imposing the 



modified field equations "phcnomcnologically" , i.e. without considering the full underlying field theory, it is 
not possible to give such a prescription. An educated guess is to impose the field equation to be causal, and 
this can be implemented by choosing the retarded propagator in the spectral representations. This is the 
correct prescription when the effective equations come from quantum corrections, and in particular for weak 
gravitational fields. In this case, this choice corresponds to the vacuum quantum state of the matter fields 
|8[ . For the rest of the paper we will follow this particular prescription, although it is clear that it is not 
the only possibility. The scale dependence of the Newton constant could be a consequence of the presence 
of extra dimensions, and in this case it is not even clear that the effective equations on the brane should 
be causal [l7j . Moreover, even if the corrections come from matter quantum fields, one could consider more 
general initial states. These issues deserve further investigation. 

Before addressing specific calculations, we introduce the following notation for the kernels: (— □ + b 2 )~ a = 



G£(x,x') and ln[ 



□ 

J? 



L(x,x'). Moreover, G(M 2 , x, x') will stand for the retarded propagator of a scalar 



field of mass M. The action of the kernels on a scalar function / is explicitly given by 



and 



d 4 x'^\^T\L(x,x')f(x') , 





□ " 




111 




'=-! 








(-□ 


+ b 2 y 





d A x'^W)\GZ(x,x')f{x') 



where 



G"(x, x') = 2sln(mr) / + °° dm m^Gim 2 + b 2 , x, x'). 
I" Jo 



(9) 



(10) 



(11) 



III. (-n)- Q IN MINKOWSKI SPACETIME 



Let us begin by computing the application of the operator defined in Eq. (|10[) in Minkowski spacetime to 
an arbitrary function. For this purpose, we start by considering the retarded propagator for a massive field, 
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which is twice the real and causal part of the Feynman propagator, 

<i 4 fc expjifc • x} 



G{m 2 + b 2 ,x,x') = 26(t-t')ft[G F (x,x')}=29(t-t')$l 



(2tt) 4 fc 2 - ie + m 2 + b 2 



(12) 



om ^ [ + °° dk y if* y/M sin(fc|f " sin(Vfc2 ± WW + fe2(t ~ *'» 

= 20(t-t) / — — ^fcexpj-e (i - i } =- ====== , 

Jo ( 27r ) \x-x'\ Vk 2 +m 2 + b 2 

where the positive number e — » at the end of the calculations. 

The generalized Green's function G% (x, x') can be then obtained by inserting this result into Eq. fTTj) and 
performing the remaining integration. In the case in which b = we get (23j 

G%(x,x') = ENpIh±2M 9 (t -H- \x-x'\)[{t t') 2 \x - x'\ 2 r~ 2 exp{- e (t - f)}, (13) 
7r z l (a — 1)1 {2a) 

which is the same as the corresponding retarded Green's function obtained by Bollini and Giambiagi in 
Ref. [24[ using a different method. Note that this generalized Green's function does not satisfy the Huygens' 
principle, i.e. its support is not restricted to the surface of the past light cone. 

With this result we can apply the (—0)~ a operator to any test function that depends on the spacctime 
coordinates. In the case of time independent functions, the time integral of the retarded propagator coincides 
with the Green's function of the Laplacian, 

+ °° A+>n< 2^,2 i\ f d 3 k expfrfc-QE-P)} 

at G(m +b ,x,x ) = / — 5 — -~ — . (14) 



(2tt) 3 k 2 + m 2 + b 2 



Therefore, in such cases we obtain 



d k exp{ik ■ x} 

-/(«)) 



(2tt) 3 (fc 2 +6 2 ) c 

where /(fc) = / d 3 a;cxp{— ifc • x}f(x). 

As one can immediately note from the equation above, the application of (— □ + b 2 )~ a to a constant 
function / = 1 gives b~ 2a . This result is the one that would be expected by expanding (— □ + b 2 )~ a in 
powers of the derivative operator □. 

On the other hand, if the test function only depends on the time coordinate we are allowed to perform 
the integration on the space coordinates, with the result 

9t{t,t)= [ d 3 x'G^x,x') =2 1 / 2 - a b 1 / 2 ~ aJ ^-J a _ 1/2 ( b ( t -t'))9(t-t')(t-t'r- 1/2 , (16) 
J T(a) 

where J is the Bessel function of the first kind [25[ . For b = it reduces to 

3o (*, = 6XP{ ~f ~ t/)} 0(t f )(* t>f«-\ (17) 
T(2a) 

As particular examples, let us consider the test functions fi(t) = (—t) 13 (t < 0,(3 < —2a) and fz(t) = 
exp{Vt} (v > 0), for which we obtain 

(-□r7i = j + ™dt>g%(t,t>)h{t>) = r( ~^~ /3) (-^) 2a+/3 1 (18) 

/+oo 
dt'gg(t,t')f 2 (t') = ^ 2Q exp{^}. (19) 
- OO 
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Note that as the definition we are considering is not applicable to all integer values of /3, one cannot obtain 
the second result from the first one by a Taylor expansion of the exponential function. Note also that both 
results can be obtained from analytic continuation, by considering the application of the derivative operator 
(— □)" with n G N that yields a real result after replacing n by —a [2(|. 

The results of this Section are useful to analyze field equations like Eq.([2]) in the weak field limit. We will 
use them in Section [VI[ to study the generation of gravitational waves in the linearized limit. 



IV. (-□)"" IN FRW SPACETIMES 

In this Section we consider three qualitatively different spatially flat FRW spacetimes. We use coordinates 
such that 



ds 2 = a 2 (r 1 ) 



= -dt 



2 , 2 

a 



(f)X>? 



(20) 



where a is the scale factor, r\ is the conformal time (drj = dt/a(t)) and t the cosmic time. 

A. De Sitter spacetime 

We set the scalar factor to be 0(77) = — (Hrf)" 1 , where the constant H = a/a (with a dot denoting 
differentiation with respect to t). The retarded propagator for massive scalar field in this spacetime is given 

by na 



G(m 2 + b 2 , x, x') = % - rfm \ i— QH, ^ F 

8tt cos{nv) 



3 3 (y-rf -ie) 2 -\x-^\ 2 

2 ^ 2 ' ' 47777' 



. (21) 



where F is the Gauss's hypergeometric function [25[ and v 2 = 9/4 — (m 2 + b 2 )/H 2 . 

As in this case the integration on the mass parameter m is very difficult, we restrict ourselves to test 
functions that depend only on the time coordinate. Then, before performing the integration on the mass we 
integrate on the space coordinates, which results in 



g(m 2 + 6 2 , v , 77') = I d 3 x'^n\G(m 2 + b 2 , x, x 1 ) = { ^ sinh f vln 1 V ' 



This equation can be integrated on the mass parameter to find 



fifa rf) = 2sin(Q ^ ) [ + °° dm m l - 2a g(m 2 + b 2 , 77, 7/). 



7T JO 

The result depends on whether v is a real or purely imaginary number. We get 



9 b_ 

4 H 



(*"») 



T(a) (-77O 5 / 2 V 2 

a N l/i-a/2 



1/4-Q/2 



5 -grin , b 2 <9/4H 



K - ^ In ( 2- 



b 2 > 9/4H 2 



(22) 



(23) 



(24) 



where / is the modified Bcsscl function of the first kind 21 
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As in Minkowski spacetime, applying the operator (— □ + b 2 )~ a to a constant / = 1 we obtain the 
"expected" result b~ 2a . As we will see in the next sections, this property is not true in spacetimes with 
initial singularities. This is important, for example, in the context of models where in order to explain the 
smallncss of the experimentally inferred dark energy density, the conventional laws of Einstein's gravity are 
modified at very large scales by introducing non local operators. In such models, gravity is modified so that 
a large vacuum energy density (identified with the dark energy density) does not give rise to large observable 
curvature. Phenomenologically, it can be seen as a result of the possibility of having a modified Einsten's 
equation like Eq.@ with the non local effective gravitational coupling such that when it is applied to a 
smooth and persistent source yields a small contribution due to G(d)l = G(0) << Gm- For example, in the 
presence of a vacuum energy density, = —p v g^, the trace of Eq.([2]) becomes 

R = 32%G(-D) Pv , (25) 

so the effective vacuum energy is expected to be G(0)p v -C GnPv As a phenomenological approach, the idea 
of addressing the cosmological problem with non local modification of gravity in the infrared is discussed for 
example in Ref. [28| , where it is argued that acausality is a fundamental feature for this program to work. 
The application of (-□)-" to f(rj) = {-Hri) k (k < 0) yields 

^'^vViv^ j^^ , (26) 

which can also be seen as an analytic continuation of the result of the application of the covariant derivative 
operator (—□)", with n an integer number, from n — > —a (we remind that we are assuming < a < I). 

It is useful to note that if we are only interested in applying the kernel (— □ + b 2 )~ a to test functions that 
are independent on the space coordinates we only need to know the zero mode of the scalar field. That is, 
we can start by working out the solution of the Klein-Gordon equation for a spatially homogeneous scalar 
field 4>(rf) = (27r) _3 / 2 x(?7)/a(r/) (with mass M = (to 2 + b) 1 / 2 and minimally coupled to gravity), 



x"(v) + a 2 (v) 



TO 2 



R 



x(v) = 0, (27) 



where a prime means a derivative with respect to r\ and R = 6a" (rj) / a 3 (rf) is the Ricci scalar. After imposing 
the Wronskian normalization condition, 

x(v)d v X* (v) - X* (v)d v x(v) = i, (28) 
the retarded propagator integrated in the space coordinates is determined by 



3 (to 2 + 6 2 ,7 7 ,t/)= / d 3 x^(^G(m 2 +6 2 ,x,xO = 20(r ] - ? /) v ^(^JR{z(2^) 3 0( ?7 )0*(7 ? ')}- (29) 
We can reobtain the result presented in Eq. (|2"2"|) by starting with the following solution to Eq. (j2"7|) : 

9[V> a (77)(27r)3/ 2 a(77)(27r)3/ 2 ] ^f., iv e M. [ ' 



As one can easily check, this solution satisfies the Wronskian condition (|28|) . Then, in order to find the 
propagator g(m 2 + b 2 ,r),r]') given in Eq. (f2"2")) we only need to insert this solution into Eq. (f2T))) . It is worth 
noting that the result does not depend on the initial condition for the mode. 

For the remaining FRW spacetimes considered below, we restrict ourselves to test functions that depend 
only on time, and therefore wc will follow this last procedure. 
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B. FRW spacetime with a(t) oc t 



1/2 



We now consider a radiation dominated FRW spacetime with scale factor a(t) = a^t 1 / 2 . This spacetime 
has an initial singularity and a finite particle horizon. 

A solution of Eq. (f2"T|) that satisfies the Wronskian normalization condition is given by 



<Kv) 



^ H[%(Mt), 



8TTa(rj) 



(31) 



where t = a^rj 2 /4, M = \/ m 2 + b 2 , and is the Hankel function of the second kind [25^ . Substituting 
this solution into Eq. (f2"5)) we obtain 



4 /7/2 

Tial r) 



g(m 2 + b 2 , V ,r,') = % - rf) ^ ' { J 1/4 (Mt) J_i /4 (Mf) - J- 1/4 (Mt) J 1/4 (Mi')} 
For 6 = the mass integration yields 

. _ _ a / „ \ 2a- 1 r 



F 



a — — , a, 2a; 1 



(32) 



(33) 



Contrarily to the previous cases, after applying the (— □) Q operator to a constant / = 1 we obtain a 
finite result, 



+°° r(5/4)< 2 «4-« 
d?7 ffo {V,V) = 



r(a + l)r(5/4 + a) 



(34) 



This finite result is a consequence of the fact that this spacetime has an initial singularity since, as it was 
pointed out in Ref.[28|], starting at t = (77 — in this case) it would take causal physics too much time to 
recognize that the constant was truly a constant . This interpretation may be further supported by computing 
the application of the (— □ + 6 2 ) _Q to/ = l, 



drf 9b(v,v') = sin(7ra)£ 2Q 2- 2Q+1 / 2 



r(-a) 



r(-l/4)r(5/4+a) 



1F2 



5 b 2 t 2 

a, 1 + a, — + a: — 

4 4 



where 1-F2 is one of the generalized hypergeometric function p F q [29(. Then, taking the limit t 
arrive at the expected result 



lim 

t — >+oo 



dv' 9"(V,v')=b- 



(35) 



-co we 



(36) 



On the other hand, despite the initial singularity, applying (— □) a to f(rj) = {a^-q 2 /4) 13 = t@ ((3 > —2) 
we get 



f +X dn' aVr, W) fU) 1-^+3 r(/3/2 + l)r(/3/2 + 5/4) 
o dV 9°W,Vmv)-4 t Tif3/2 + a + i)r ((3/2 + 5/4 - 



+ a) 



(37) 



This result can also be obtained through analytic continuation of the application of covariant derivative 
operator (— for positive integers n [15| . 
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C. FRW spacetime with a(t) cc t 2 



We will now consider a FRW spacetime with an initial singularity but. in contrast to the previous example, 
without a finite particle horizon. The choice a(t) = a^t 2 will allow us to find a solution of Eq. (j2"7| in terms 
of elementary functions. Indeed, a solution that satisfies the Wronskian condition in Eq. (j2"5)) is given by 



exp{— iMt} 



M 



a( V )(2ir) 3 / 2 (Mt) 3 V 2a 



\MH 2 - 3(1+ iMt)], 



with t = —(aQi]) 1 and M = V rri 2 + b 2 . Inserting this solution into Eq. (|2"9"l) we find 



g{m 2 + b 2 , 77,77') = 



jcos(M(f - *')) 3Mt(M 2 t' 2 - 3) - 3Mt'(M 2 t 2 - 3) 



a MHH' 3 a(rj) 
+ sin(M(t-f')) 9M 2 tt' + (M 2 t 2 - 3)(M 2 t' 2 - 3)J | , 

and for 6 = 0, after integrating on the mass parameter, we obtain 



(38) 



(39) 



9o(V,v') 



a( v ')6( V -r)')(t-t') 



l\2a-l 



t' t' 

3- + 3— + 3(2a- 1) 
t t 5 v 



T(2a)(l + 2a)(3 + 2a) 
The application to a constant / = 1 yields the result 

H 9o(v,v') 



'*1 
t 2 + 1* 



+ (3 + 4a(o-l)) 



15* 



2n 



T(2a + 2)(2a + 3)(2a + 5)' 

which is also finite. 

As in the previous case, applying (— □ + b 2 )~ a to a constant we get 

b 2 t 2 ' 



15 I 

drf 9b(v,v') = — * 2a sin(27ra) j T(-3 - 2a)iF 2 

+ 3r(-4-2a)iF 2 
and taking the limit * — > +00 we find 



a, 2 + a„ - + a: 



(40) 



(41) 



(42) 



5 6 2 * 2 
a, - +a,3 + a; — 



3r(-5 - 2a)!F 2 



b 2 t 2 



a, 3 + a, 



lim 

£ — >+oo 



(43) 



Therefore, the fact that the application of the (— D)~ a to a constant yields a finite result seems not related 
to the existence of a particle horizon, but to the existence of a singularity at * = 0. 

Another interesting property we can point out is that considering the advanced propagator, 



gadv(m 2 + b 2 ,7 h 7/) = -20(7/ - 77)^7^3? {7(27r) 3 0(77)0* (77')} , 

instead of the retarded one, we obtain 

(-v + b 2 u v i = b~ 2a , 



(44) 



(45) 



independently of the value of the initial time t. 

Finally, for the sake of completeness, we compute the application of (— D) - " to f(r]) = (—aQi])^ 13 = t' 3 
(f3 > -2), with the result 



f\ _ A -a.2a+0 



r(/3/2 + 7/2)r(l + /3/2) 



r(l + a + /?/2)r(a + 7/2 + 0/2) : 
which can also be obtained by analytic continuation from an integer number n to —a. 



(46) 
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V. THE OPERATOR ln[ D/V] 



In this Section we compute the application of the operator that can represent the logarithm of the covariant 
d'Alembertian to scalar test functions that are independent of the space coordinates, for each of the four 
spacctimcs considered above. As described in Section UIl we will do that by taking the limit a — > in the 
kernel (-□)"". 

We begin with Minkowski spacetime. We start by applying the kernel (— n)~ Q , given in Eq. (fT"7[) , to a 
generic test function / that is independent on the space coordinates, 



+00 



dt' v 2a g%(t,t')f(t') 



r(2a) 



dt' n 2a cxp{-e(t - t')}(t - t') /(O 



(47) 



where we have added a constant fi with mass dimensionality. 

Before expanding in powers of a, and in order to exclude explicitly the local contribution, we perform an 
integration by parts, 



1 



r(2a) 



/ ..2 Q exp{-e(i-i')} 



dt' [i 



/(*') = fit) 



1 



W) J_^^T(2a,e { t-t'))- m , (48) 



where T(a, x) is the incomplete gamma function [25[, and we have assumed that f(t')T(2a, e(t — t')) goes to 
zero as t' — ► —00. 

After subtracting the local zeroth order contribution f(t) and expanding up to first order in a we find 



+00 



dt' gg(tj)f(t') - f(t) ~ -ai\n(j^j + 2 J* dt'T(0,e(t - t'))-^f(t') 



Finally, in the limit e^flwe obtain 



In 



□ 



I = -2 7 /(i) -2 I dt' In (n(t - t'j) — f(t'), 



(49) 



(50) 



where 7 is the Euler's constant and we have assumed that f(t') — > as t' — > —00. The same result can also 
be obtained from the representation proposed for this kernel in Ref. [30l | . 

With the use of this integral representation one can straightforwardly compute the application of ln[— 

to spatially homogeneous functions. For example, for f(t) = (—t)" (with t < and j3 < —2a) it yields 



In 



□ 



f=(-tf{-\n(t 2 ^) + 2^(-(3)}, 



(51) 



where ip(x) is the Psi (Digamma) function [25|]. Note that it can also be obtained by expanding the result 
in Eq. psp in powers of a and extracting the first order contribution. 

For the three remaining spacetimes, we can follow the same procedure. Details of the calculations are 
relegated to the Appendix lAl 

For the dc Sitter spacetime, we arrive at 



In 



□ 



/ = In 



H- 



f(rj) - drf In 



In 



Ei -3 In 



dr) 



7 fW), (52) 



where Ei is the exponential-integral function [231 ] and we have assumed that f(j]) — > as ?y — > —00. The 
application of this operator to f(rf) = (—Hrj) k (k < 0) gives the result 



In 



□ 



(53) 
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Likewise, for the radiation dominated FRW spacetime (a(t) cc i 1 / 2 ) we obtain 



In 



□ " 






f = 


7V\ 





4-27 In 

' 2 



In 1 



drf 



4rf- 
5?7 5 



F 



'4' 4' t? 4 



(54) 




with /(?/) — > as rj — > 0. After applying it to / = t' 3 = {a -q/2) 213 (j3 > 0) we get 



In 



/' z 



/ = if In 



16i2"' 



( , 



/9 



V' TT + 7 + 77 + 1 



(55) 



Lastly, in the case of the FRW spacetime with a(t) oc t 2 the result is 



In 



/ = 



4G 



2t'° 



2f_ 

T 



■it' 



m-l dt> t iW + - + — + 2ln(l- T )}-f(t>), (56) 



3t 3 




with f(t) -> as t 0. For / = ^ = (-ao??) _/3 (/? > 0) it yields 



In 



□ 

T 2 



f = t \ In 



m 2 



+ v 







1 



(57) 



As we have pointed out in the Introduction, there are different proposals in the literature to take into 
account the running of coupling constants, which correspond to different replacements of the energy scale ft 
by 1/t, H, p,/a(t), etc. In cosmological situations, it is common to assume a power law behavior for both the 
scale factor and the matter content. In such situations, the relevant results arc the ones given in Eqs. (f57|) . 
(|4"6")h (|5"5")) . and ([ST)) , which were obtained by applying the corresponding non local kernel to test functions 
that have a power law dependence on the time coordinate (cither r\ or t). These results can be written in 
the form 



(-□)-°7 ex H- 2a f 



In 



□ 

7 



[ ^2 ) 



(58) 
(59) 



Moreover, in de Sitter spacetime, if we consider test functions that have a power law dependence on the 
conformal time 77, we can also write the results in the same form (sec Eqs. ([26|) and l|53[)). This would suggest 
that in FRW spacetimes the action of the non local kernel is equivalent to the replacement p, — ► H. However, 
this is not strictly correct. Indeed, the results depend on the test function. On the one hand, the constant 
of proportionality in Eq. (|58[) depends on a and (3. On the other hand, in de Sitter spacetime, if we apply 
the non local operator to a function / = (— t)* 3 (J3 < and t < 0) we obtain 



In 



□ 



(60) 



where F(Ht,(3) is a complicated function whose functional dependence on time is not logarithmic. For 
example, in the particular case of (3 = —2 it reduces to 



F(Ht, -2) = 2 - 7 - M-Ht/3) - exp{-3Ht}(l + 3Ht)Ei(3Ht). 



(61) 



This non logarithmic behavior may have been expected from the fact that, in contrast to the other FRW 
examples, for de Sitter spacetime one has t and H as two different and relevant scales. 
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VI. GENERATION OF GRAVITATIONAL WAVES 



In this Section wc analyze the modified Einstein's equations (|2|) with a running coupling given by Eq. (J6)), 
for the particular case b = 0. We assume weak gravitational fields. Using an expansion in powers of the 
curvature, it can be shown that the modified equations follow from the non local action [13] 



S = ~J <r X y/\g\ {Jq^ - ) D R ^ + °i R U)) + Matter, (62) 

where S ma tter corresponds to the source of the gravitational waves {25S ma tter = \/W\T tlJ ij8g^ y ). Here, the 
kernel F[— □] = (— L 2 n)~ a represents a small modification of Einstein's theory for very large scales. 

As anticipated, we work in the weak field limit Q^i* - — ?7/iv + h^, where describes a linear perturbation 
to the Minkowski metric. In the harmonic gauge, up to first order in h^, the modified Einstein's equations 
become 

(l + F[-n])nh lu , = -l&rG N T MV , (63) 

where = h^, - 

We assume that the inclusion of the kernel F [— □] yields a perturbatively small correction to the standard 
result. Therefore, we split = h ®} + h^J, where h!^} satisfy 

nh$ = -WirG N T^, (64) 
and hflu can be perturbatively computed by solving the following equation: 

D ^ = -(-L 2 n)- a ah$ = 16TrG N (-L 2 D)- a T^ = -16ttG n T^. (65) 
With the use of the retarded Green's function we can write 

h$(t,x) = 4G N [ -A^f^t - \x-x'\,x'), (66) 
J \x — x'\ 

where 

f MV (t,sS) = - J d 4 x'G%(x,x')T^(x'). (67) 

Note that the value of hfy in the spacetime point x is determined only by the values that takes on the 
past light cone of x, but in order to compute one needs to know all the values of the energy-momentum 
tensor inside that past light cone. Therefore, the solution for the metric perturbation does not obey 
the Huygens' principle. This last issue was analyzed in detail in Ref.[24| for different spacetime dimensions. 

Finally, working with the Fourier transform and using the generalized Green's function of Eq. (|13|) one can 
find 

k) = -i6nc N {i + J£ ) 2 + ifci 2 rn - T ^: k \ rm - m 

[-(uj + ie) 2 + \k\ 2 ] 

The non locality of the modified Einstein equations has other interesting consequences. It gives rise to 
corrections to the Newtonian potential which in turn imply that the metric outside a massive object may 
depend on its inner structure, violating Birkhoff theorem [31(. Therefore, we expect the corrections to the 
metric outside a collapsing star to be time dependent, even in spherically symmetric situations. 
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VII. CONCLUSIONS 

We have presented a detailed derivation of the action of the non local kernels (— D)~ a and ln[— □] on scalar 
functions. Our starting point was the representation of both kernels in terms of integrals of the retarded 
propagator of a massive scalar field. 

In Minkowski spacetime, we derived an explicit expression for the kernel (— D)~ a . This is relevant for 
the analysis of the modified Einstein Eq.([2]) in the weak field limit. As an application, we obtained formal 
expressions for the generation of gravitational waves, and showed explicitly that Huygcns' principle is not 
longer valid when one considers a scale dependence of Newton constant. 

In Robertson Walker spacetimes, wc considered the action of (— D) - " on spatially homogeneous time 
dependent test functions. We presented a method for obtaining general expressions for the action of the non 
local kernel on such test functions, which only involves the zero mode solution of the Klein-Gordon equation 
for a massive, minimally coupled scalar field. We specialized the general expressions to particular cases, as 
test functions defined as powers of the cosmic or conformal times. We have also obtained an expression for 
the action of ln[— □] on time dependent test functions by considering the limit a — > 0. 

The application of the operator (— D)~ a to functions that depend on t in Robertson Walker spacetimes, 
or functions that depend on r in static metrics with spherical symmetry, has been considered previously by 
Hamber and Williams in Refs.pjl [2lJ. Their approach was based on an analytic continuation of the action 
of the derivative operator (— □)" for positive integers n to the case n — > —a. For simple test functions in 
which the analytic continuation can be carried out, this procedure gives the same answers than the non 
local kernels defined in this paper. However, it is worth to note that, within our approach based on the 
integral representation (Eq. <J7J>), it is in principle possible to evaluate (—□)""/ for arbitrary scalar functions 
in arbitrary spacetimes. Moreover, this representation shows explicitly the non local properties of the kernel 
and could also be generalized to tensor functions. 

In the cases of the FRW spacetimes in which both the scale factor and the test function have a power law 
dependence on either the cosmic or conformal time, after applying the non local kernels, we have found that 
the functional dependence on time of the results is the same as the one obtained after the replacement of 
the energy scale \x by the Hubble rate H. However, using as an example the de Sitter spacetime, we have 
explicitly shown that in general the result depend on the test function and can be very different from the 
one obtained with the use of any of the replacements mentioned in the Introduction. 

We also investigated whether the action of 1/G(— □) on a constant test function gives 1/G(0) or not, 
as naively expected for analytic functions of the d'Alcmbcrtian. In particular, we have shown that ( — □ + 
b 2 )~ a l = b~ 2a only for spacetimes without initial singularities, or when considering acausal propagators, 
and that it seems not to be related to the existence or not of a particle horizon. These results are in tune 
with the claim p8j that, in the presence of an initial singularity, acausality may be crucial to solve the 
cosmological constant problem using infrared modifications to General Relativity. 

As for future work, we consider that a similar method as the one proposed in this paper could also be 
used, without much additional difficulty, in static and spherically symmetric spacetimes for test functions 
that respect such symmetries. 
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APPENDIX A: EXPANSION UP TO FIRST ORDER IN a 

In this appendix we describe some technical details about the computation of the expansion up to first 
order in a of the application of (— D)~ a to time dependent functions. The procedure we follow here is the 
same as the one described in the main text for the Minkowski case. 

1. De Sitter spacetime 

In order to obtain the result presented in Eq. (|52]) . we start from Eg. ([24]) (in which the zero mode Green's 
function g^ is given). We compute a primitive of the Green's function g$ with respect to rf , 

r( i\ f,'2 Qa( i \ 3 V 2 -V« r(l/2) n' drf' / ?/ \ 3 / 2 ^-i/ 2 (|ln(f)) 
Q(V,V) = / drf n 9 (V,V) = 77^ WZT / — 177 37 ; ; 7TT/2^ ( A1 ) 



H 2a T(a) J -rj" \rf' J (, fr,"\\ 1/2 ~ a 



3 ^-v« r(i/ 2 , h\ w U±\y> n {^ w (* u (±\\ +I „ +1 Jiut 



2(H) 2a T(a + l)\n'J V \V J J { a ~ 1,z \2 \r)JJ a ^^\2 \r] 

where the second equality follows after the change of variable s = 3/2 \\\{rf' /rf) (for the last equality see 
Ref.0). 

Assuming that G(j], rf)f(rf) ~ (}n{rf / rf)) a f (jf) goes to zero as rf — > —00 (a — > + ), and expanding (77, rj') 
up to first order in a, 

' e 1 p? f ( ( rf 

y{'h >l ) - -J- - « "j I 

we get 



0fo )f /)~-l-a|]n(^]n(i] ) + ^ ( -31n ( ^ ) ]>, 



<V M aa S?fa, »/)/(»/) - /fa) ^ « */ {in (|^£ln + (-3 In (A3) 

where is the exponential-integral function (231 ] . Then, the required operator (|52p is obtaining by comparing 
this result with the first order of the alpha expansion of (— D)~ a . 

2. FRW spacetime with a(t) oc t 1/2 
We begin by computing a primitive of the Green's function g$ given in Eq. (|55|) (see Ref. [HI), 

,4" 

(A4) 



qm) - /^^v)=4-^/;'^(i-^i f 



2a-l r a 

a — — , a, 2a; 1 j- 



2n 



r(5/4)r( a + i/2)(^) 2 " (^) 2 « ^_!/!| 



r(a + 5/4)r(l/2)r(2a+l) 4«r(2a + l) \ r/ 4 
Then, the integration by part results in 



, , , - a 

4 rf 



1 r/ 4 

a — -, a, 2a + 1; 1 

4 ?y 4 



dv'f^ a 9S(v,rf)f(rf) = GMfiv) - drfg{r,,ff)-^f{ff), (A5) 
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where we have assumed that f(rj')0(rj,r]') ~ f(r]')ri' — > as ?/ — > 0. It is straightforward to calculate the 
expansion up to first order in a of Q{rj, rj), 



Q(t], rj) a 1 + a(-4 + 2 7 + - + ln(2/i 2 t 2 )). 
In order to expand Q{rj, rf) it is useful to use the following identities (see Ref.[25[): 



F 



a — -, a, 2a + 1; 1 — z 
4 



r(2a + l) f T(5/4) 



r(a+l) lr(a + 5/4) 

5/4 r(-5/4) f 



F 



1 1 



a z 



F 



1 1 

a , a, ; z 

4' ' 4' 



r(a- 1/4) 
1 + a(l - 4a) / ds F 



5 , 1 
a + -, a + 1, -;z 
4 4 



3 , 3 
Q'H — , a + 1 , — : s 

4' '4' 



Jo 



1 + a ds F 



3 3 

— [ — • e 

4' '4' 



= 1 - aln(l - z). 



Then, we find that up to first order in a 



(A6) 



(A7) 



(A8) 



5(77,77') ~ a 



5?7 5 



5 9 rf_ 

4' '4 ; ry 4 



In 1 



(A9) 



Inserting (|A6|) and (|A9|) into Eq. (|A5[) . and comparing with the expansion of (— □) a we arrive at the result 
given in Eg. ([54]) . 



3. FRW spacetime with a(t) a t 2 
A primitive of the Green's function given in Eq. (|40p can be obtained as 

(t-t')2«-l M 2 Q 



g(r,, V ') = / drf n 2a gS( V , V ') 



t' J 



dt r(2a)(l + 2a)(3 + 2a) 
't' 2 t' 4X 



(A10) 



3- + 3-=- + 3(2a - 1) — + — + (3 + 4a(a - 1)) 



t t 



t 2 t 4 



t 3 



where we have used that a(rf)drf = dt' . In this case, it is not difficult to carry out the a expansion to find 
that up to first order we have 



-^+ln( e 2 V^ 2 ) 



5(77,77') ~ a { — - - 



2f 2t' 5 2t ,s 



t 5t 5 3t 3 

where we have assumed that G(r),r)')f(tf) ~ t' 2 f(t') -> as i! -> 



- 2 In 1 - 



(AH) 
(A12) 
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